RAY  W-  HERRICK  LABORATORIES 

A  Graduate  Research  Facility 
of  The  School  of  Mechanical  Engineering 


PERFORMANCE  OF  STRAIGHT-FIN  AND 
MICROCHANNEL  HEAT  EXCHANGERS  IN 
STEADY  AND  PERIODIC  FLOWS 


Sponsored  by 
Office  of  Naval  Research 


HL  2001-18  Report  #1614-1 


Purdue  University 


West  Lafayette,  Indiana  47907-1077 


PERFORMANCE  OF  STRAIGHT-FIN  AND 
MICROCHANNEL  HEAT  EXCHANGERS  IN 
STEADY  AND  PERIODIC  FLOWS 


Submitted  by: 


Approved  by; 


Sponsored  by 
Office  of  Naval  Research 


HL  2001-18  Report  #1614-1 


Aaron  Alexander,  Graduate  Research  Assistant 
Luc  Mongeau,  Co-Principal  Investigator 
James  Braun,  Co-Principal  Investigator 

Robert  J.  Bernhard,  Director 
Ray  W.  Herrick  Laboratories 


OCTOBER  2001 


20011018  038 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  lor  this  collection  of  information  is  estinuted  to  average  1  hour  per  response,  including  the  time  for  reviewing  insUuclions.  searching  data  sources, 
gathering  and  maintaining  the  data  needed,  and  complclcng  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Service.  Directorate  for  Information  Operations  and  Reports. 

1215  Jefferson  Oavis  Higfr^y,  Suite  1204.  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget. 

Paperwork  Reduction  Pi^ect  (0704-0188)  Washington,  DC  20503. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

08-10-2001 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  DATE 

Final  Report 


3.  DATES  COVERED  (From  -  To) 

May  1999  -  May  2001 


5a.  CONTRACT  NUMBER 


PERFORMANCE  OF  STRAIGHT-FIN  AND  MICROCHANNEL 
HEAT  EXCHANGERS  IN  STEADY  AND  PERIODIC  FLOWS 


6.  AUTHOR(S) 


Alexander,  Aaron 
Mongeau,  Luc 
Braun,  James  E. 


5b.  GRANT  NUMBER 

N00014-99-10642 


5c.  PROGRAM  ELEMENT  NUMBER 

PE  61153N 


5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
Purdue  Research  Foundation 
Division  of  Sponsored  Programs 
Hovde  Hall,  Third  Floor 
West  Lafayette,  IN  47907 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 
Office  of  Naval  Research,  ONT  331 
Ballston  Centre  Tower  One 
800  North  Quincy  Street 
Arlington,  VA  22217-5600 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

HL2001-18 


110.  SPOUSOR/MONiTOR’S  ACRONYM(S) 


11.  SPONSORING/MONITORING 
IGEWCY  REPORT  NUMBER 


12.  DISTRIBUTION  AVAILABILITY  STATEMENT 


Approved  for  public  release:  distribution  unlimited 


14.  ABSTRACT 

The  performance  of  straight-fm  and  microchannel  heat  exchangers  was  investigated  in  steady  and  periodic 
flows.  Steady  flow  measurements  were  performed  using  an  open  test  section  wind  tunnel.  Periodic  flow 
measurements  were  performed  using  a  thermoacoustic  cooler  prototype.  Detailed  temperature  measurements  were 
made  for  varying  operating  conditions  using  thermocouple  arrays.  The  results  suggest  that  nondimensional  heat 
transfer  coefficients  based  on  stack  end  temperatures  calculated  from  steady  flow  data  may  be  applicable  to  predict 
the  thermal  performance  in  a  periodic  flow,  provided  the  root-mean-square  local  particle  velocity  is  used  as  the 
velocity  scale.  The  accuracy  of  rms-based  steady  flow  correlations  may  be  comparable  or  in  some  cases  better  than 
that  of  linearized  boundary  layer  models  frequently  used  in  thermoacoustics. 


15.  SUBJECT  TERMS 

Heat  exchangers,  thermoacoustics 


16.  SECURITY  CLASSIFICATION  OF: 


b.  ABSTRACT 

c.  THIS  PAGE 

u 

u 

17.  LIMITATION  OF 
ABSTRACT 


18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
OF  PAGES 


19b.  TELEPONE  NUMBER  [Include  ana  code) 


Standard  Form  298  (Rev.  8-98) 
Prescribed  byANSI-SldZ39-18 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES . 5 

NOMENCLATURE . 6 

1.0  INTRODUCTION  11 

2.0  HEAT  EXCHANGER  DESCRIPTION . 13 

2.1  Straight  Plate  Finned-Tube  Heat  Exchanger . 13 

2.2  MicroChannel  Heat  Exchanger . 15 

3.0  EXPERIMENTAL  APPARATUS . 16 

3.1  Wind  Tunnel  for  steady  flow  measurements . 16 

3.2  Thermoacoustic  System  for  Unsteady  Flow  Measurements . 18 

4.0  EXPERIMENTAL  PROCEDURES . 20 

4. 1  Steady  Flow  Experimental  Procedures . 20 

4.2  Oscillating  Flow  Experimental  Procedures . 20 

5.0  ESTIMATION  OF  HEAT  TRANSFER  COEFFICIENTS . 21 

5.1  Steady  Flow  Calculations . 22 

5.1.1  Theoretical  Calculations . 22 

5.1.2  Calculations  from  data . 23 

5.2  Description  of  Oscillating  Flow  Calculations . 27 

5.2.1  Theoretical  Calculations . 27 

5.2.2  Calculations  from  Data . 29 

6.0  RESULTS . 30 

6.1  Steady  Flow  Results . 30 

6.2  Performance  of  the  Heat  Exchangers  in  Oscillating  Flow . 32 


3 


7.0  CONCLUSIONS  AND  FUTURE  WORK 


38 


LIST  OF  REFERENCES . 39 

APPENDICES . 41 

I:  Table  of  Thermoacoustic  Raw  Data  for  the  MicroChannel  Heat  Exchangers . 41 

II:  Table  of  Thermoacoustic  Results  for  the  MicroChannel  Heat  Exchangers . 44 

III:  Table  of  Thermoacoustic  He-Ar  Properties . 46 

IV:  Table  of  Water-Side  Properties . 48 

V:  Table  of  Thermoacoustic  Raw  Data  for  the  MicroChannel  Heat  Exchangers . 50 

VI:  Table  of  Thermoacoustic  Results  for  the  MicroChannel  Heat  Exchangers . 50 

VII:  Table  of  Thermoacoustic  He-Ar  Properties . 51 

VIII:  Table  of  Water-Side  Properties . 5 1 

IX:  Thermoacoustic  Sample  Calculation . 52 


4 


LIST  OF  FIGURES 


Figure  Page 

1  Drawing  of  Fin-tube  Heat  Exchanger  14 

2  Picture  of  Installed  Fin-Tube  Heat  Exchanger  14 

3  Drawing  of  MicroChannel  Heat  Exchanger  15 

4  Picture  of  Installed  MicroChannel  Heat  Exchanger  16 

5  Schematic  of  the  Wind  Tunnel  used  for  the  Steady  Flow  Measurements  17 

6  Picture  of  the  Wind  Tunnel  used  for  the  steady  Flow  Measurements  17 

7  Schematic  of  Thermoacoustic  Cooler  19 

8  Picture  of  the  Thermoacoustic  Cooler  19 

9  Steady  Flow  j-Colbum  vs.  Reynolds  Number  31 

10  Steady  Flow  Heat  Transfer  Coefficient  vs.  Reynolds  Number  31 

1 1  J-Colbum  vs.  Reynolds  Number  for  Periodic  Flow;  32 

Fin-Tube  Heat  Exchanger 

12  Heat  Transfer  Coefficient  vs.  Reynolds  Number  for  Periodic  Flow;  33 

Fin-Tube  Heat  Exchanger 

13  J-Colbum  vs.  Reynolds  Number  for  Periodic  Flow;  35 

MicroChannel 

14  Heat  Transfer  Coefficient  vs.  Reynolds  Number  for  Periodic  Flow;  36 

MicroChannel 

15  J-Colbum  vs.  Reynolds  Number;  MicroChannel  Heat  Exchanger;  36 

16  j-Colbum  vs.  Acoustic  Reynolds  Number;  MicroChannel  Heat  Exchanger  37 

17  j-Colbum  vs.  Acoustic  Reynolds  Number:  Comparison  of  Heat  Exchangers  37 


5 


NOMENCLATURE 


English  Symbols 

Af  =  entire  fin  surface  area,  (m^) 

A  =  total  heat  exchanger  surface  area,  A  =  A^+  Aj-  (m^) 

Afc  =  fin  cross-sectional  area,  =  wt ,  (m^) 

Ao  =  external  total  surface  area,  (m^) 

Ac  =  flow  area(air  side  flow),  (m^) 

A„=  tube  internal  cross  sectional  area  (m^) 

A,  ~  entire  tube  external  surface  area  (m^) 

A/  =  tube  internal  surface  area,  Aj  ,  (m^) 

Cp=  constant  pressure  specific  heat,  (J/kg-K) 

4LA 

Dh  =  hydraulic  diameter  (air-side  flow),  = - - ,  (m) 

A 

Di  =  secondary  flow  tube  inside  diameter,  (m) 

Do  =  secondary  flow  tube  outside  diameter,  (m) 

/g  =  friction  factor  for  flow  in  smooth  tubes,  f^  =  (.791nReD  -1.64)“^ 

ID  A 

Gz=  Graetz  number,  Gz  =  — 2 — ^ 

(L/Di) 

h  =  heat  transfer  coefficient,  (W/m^'K) 

hi  =  internal  flow  heat  transfer  coefficient,  (W/m^'K) 

ho  =  external  flow  heat  transfer  coefficient,  (W/m^  K) 

j  =  j-Colbum  factor,  dimensionless  heat  transfer  coefficient,  j  =  St  Pr 
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k  =  thermal  conductivity,  (W/m-K) 

ksf^  thermal  conductivity  of  secondary  fluid,  (W/m-K) 

L=  fin  length,  length  of  flow  path  of  interest,  (m) 

LMTD  =  log  mean  temperature  difference,  (K) 

L,=total  tube  length,  (m) 

m  =  mass  flow  rate,  (kg/s) 


2  hP  .  .2, 
m  = - ,(m  ) 


kA, 


A/t=number  of  tubes  in  exchanger 


hD 

Nu  =  Nusselt  number,  Nu  =  — 

k 


P  =  perimeter,  P  =  2w  -t-  2t ,  (m) 


c  p 

Pr  =  Prandtl  number,  Pr  =  — ^ 

k 

Q  =  heat  transfer  rate,  (W) 

R  =  radius  of  heat  exchanger  viewed  by  flow,  (m) 

Rc,i=  secondary  tube  internal  convective  resistance,  (K/W) 

Re  =  Reynolds  number.  Re  = 

P 

Rfins  =  combined  fin  conductive  and  convective  resistance,  (K/W) 

Rpj  =  convective  resistance  between  secondary  tube  wall  and  primary  fluid,  (K/W) 
Rtw  =  secondary  tube  wall  conductive  resistance,  (K/W) 

Si  =  tube  spacing  in  the  air  flow  direction,  (m) 

Nu 

St  =  Stanton  number,  St  = - 

RePr 
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St  =  tube  spacing  normal  to  the  flow,  (m) 
s  =  spacing  between  adjacent  fins,  (m) 

T  =  temperature,  (K) 

r5=distance  between  tube  centerlines,  (m) 
t  =  fin  thickness,  (m) 

tt=  tube  height  (microchannel  heat  exchanger),  (m) 

UA  =  overall  system  heat  transfer  coefficient,  (W/K) 
w  =  fin  depth  in  direction  of  flow,  (m) 

Wt=  tube  width  (microchannel  heat  exchanger),  (m) 

Greek  Symbols 

a  =  thermal  diffusivity  =  — ^ ,  (mVs) 

pCp 

St  =  thermal  penetration  depth,  (m) 

rjo  =  overall  fin  efficiency,  ri^  =  1 - ^(1  -  tif ) 

A 

T]f  =  single  fin  efficiency,  rj.  (for  straight  or  pin  fin) 

mL 

//=  dynamic  viscosity,  (N  -si m^} 

(0  =  frequency,  Hz. 

Subcripts 

h  =  hot 

c  =  cold,  convective 
i  =  in,  internal 
o  =  out,  outer 
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t  =  tube 


tw  =  tube  wall 
sf  =  secondary  fluid 
5  =  surface 
fc  =  fin  cross-section 


1.0  Introduction 


Thermoacoustic  cooling  involves  heat  transfer  between  gas  particles  within  high 
amplitude  acoustic  standing  waves  and  heat  exchangers  immersed  into  the  working  gas.  The 
dynamic  pressure  and  velocity  changes  inherent  in  an  acoustic  wave  cause  a  net  flow  of  heat  up  a 
mean  temperature  gradient  along  the  direction  of  sound  wave  propagation. 

The  so-called  thermoacoustic  refrigeration  cycle  describes  the  heat  transfer  process 
within  the  working  gas,  which  is  comparable  to  a  reverse  Brayton  cycle.  Thermoacoustic  heat 
pumping  systems  include  a  porous  substrate  called  a  stack,  which  enhances  the  heat  flux  within 
the  working  gas.  Heat  is  rejected  to  and  pumped  from  hot  and  cold  heat  exchangers,  which  are 
usually  fin-tube  and  are  generally  placed  on  each  end  of  the  stack.  The  cold  heat  exchanger 
provides  heat  to  the  working  gas  removing  it  from  the  secondary  fluid  that  flows  through  the  heat 
exchanger  channels.  Conversely,  the  hot  heat  exchanger  removes  heat  from  the  working  gas  and 
rejects  it  into  another  secondary  fluid  loop  (Swift,  1999). 

Thermoacoustic  particle  displacement  amplitudes  are  typically  between  1  and  10  mm, 
and  the  operating  temperature  differences  between  the  working  gas  and  the  secondary  fluid  are 
relatively  small  (<  5  °C).  Proper  design  and  sizing  of  the  heat  exchanger  is  of  paramount 
importance.  This  includes  proper  fm  spacing,  tube  design,  tube  spacing,  and  fin  stream  wise 
depth. 

In  order  to  predict  the  performance  of  heat  exchangers  at  the  design  stage,  before 
implementation  in  a  working  system,  analytical  or  semi -empirical  models  are  necessary.  These 
models  must  account  for  the  fact  that  there  is  no  forced  mean  flow,  and  oscillatory  gas  pulsations 
are  present  within  the  thermoacoustic  environment.  It  is  necessary  to  estimate  heat  transfer 
coefficients  for  oscillating  flow  conditions. 
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Most  heat  exchanger  models  reported  in  the  literature  assume  steady  flow  conditions  in 
the  primary  working  fluid.  Currently,  there  are  very  few  tools  for  predicting  heat  transfer  in 
oscillating  flows.  Mozurkewich  (1995,1998a,  1998b,  2000)  has  developed  numerical  prediction 
tools  for  thermoacoustic  design.  Additionally,  he  has  conducted  experiments  on  thermoacoustic 
heat  transfer  from  a  cylinder,  and  from  transverse  tubes  in  pulsatile  flow.  He  found  that  for  large 
acoustic  velocity,  the  heat  transfer  in  acoustically  driven  flow  can  be  approximated  by  a  time- 
averaged  laminar  steady  flow  correlation.  Swift  (1992,  1999)  has  also  conducted  many  analytical 
and  experimental  studies  of  thermoacoustic  phoenomena,  including  thermoacoustic  heat  transfer 
this  general  area  in  which  he  was  also  looking  to  better  predict  thermoacoustic  heat  transfer  for 
design.  In  his  thermoacoustic  handbook  (1999),  the  general  ideas  for  various  heat  exchangers 
used  in  thermoacoustic  devices  are  given.  There  have  been  a  few  other  noteworthy  contributions 
to  this  field,  including  Wetzel  (1999,  2000)  who  studied  the  effect  of  heat  transfer  from  a  single 
plate  and  Brewster  (1997)  who  studied  the  effects  of  heat  transfer  between  the  elements  in  a 
thermoacoustic  system.  Yet,  very  few  experimental  studies  for  actual  heat  exchangers  in-situ 
have  been  conducted.  Mozurkewich’s  experiment  was  done  using  a  wire  in  a  tube  and  Wetzel’s 
was  done  with  a  single  plate.  Their  results  cannot  be  directly  applied  to  find  a  correlation  in 
oscillating  flow  in  real  thermoacoustic  systems  due  to  the  different  geometry.  Therefore, 
experiments  for  real  heat  exchangers  are  valueable  and  should  be  done.  The  objectives  of  the 
work  presented  in  this  report  were  to  develop  methods  for  predicting  the  performance  of  heat 
exchangers  in  an  oscillating  flow  environment,  and  to  investigate  the  performance  of 
microchannel  heat  exchangers.  In  order  to  accomplish  these  goals,  both  steady  and  oscillating 
flow  tests  were  performed  for  one  conventional  fin-tube  and  one  microchannel  heat  exchanger. 
The  fin-tube  heat  exchanger  was  designed,  built,  and  installed  in  a  thermoacoustic  prototype 
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designed  for  a  low-capacity  air-conditioning  application  in  a  previous  study  (Minner  et  al., 

2000).  The  microchannel  heat  exchanger  was  designed  and  built  specifically  for  the  purpose  of 
this  study.  Both  heat  exchangers  were  installed  in  similar  flanges,  which  facilitated  their 
installation  in  the  thermoacoustic  system.  The  existing  thermoacoustic  prototype  was  used  as  the 
environment  for  testing  the  heat  exchangers  in  an  oscillating  flow.  A  small-scale  wind  tunnel 
was  also  developed  for  testing  the  steady  flow  performance  of  both  heat  exchangers. 

There  were  several  reasons  for  performing  steady  flow  tests,  including:  1)  establishing 
procedures  for  estimating  heat  transfer  coefficients  from  measurements  and  for  correlating  heat 
transfer  coefficients  with  dimensionless  parameters,  2)  providing  baseline  data  for  comparing  the 
performance  of  the  different  heat  exchangers  under  identical  flow  conditions,  and  3)  validating 
theoretical  steady  flow  heat  transfer  coefficient  correlations  that  could  be  adapted  for  use  in 
estimating  heat  transfer  coefficients  for  oscillating  flow  environments. 

For  oscillating  flows,  special  procedures  were  established  for  estimating  heat  transfer 
coefficients  from  measured  data  and  for  correlating  the  heat  transfer  coefficients  with 
dimensionless  parameters.  Heat  exchanger  specific  correlations  are  appropriate  for  use  in  system 
simulations.  However,  they  are  not  useful  in  predicting  the  effects  of  heat  exchanger  design 
changes.  For  this  purpose,  two  different  approaches  were  investigated  for  estimating  heat 
transfer  coefficients  using  design  parameters.  One  approach  involves  a  theoretical  calculation  of 
conduction  across  a  laminar  boundary  layer  that  is  used  in  the  DeltaE  computer  program  (Ward 
1993).  The  second  approach  involves  adapting  existing  steady  flow  heat  transfer  correlations  to 
oscillating  flow  environments. 
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2.0  Heat  Exchanger  Descriptions 

Two  different  heat  exchangers  were  investigated  in  this  study.  Firstly,  a  conventional  fin- 
tube  heat  exchanger  was  used  for  preliminary  studies  and  to  develop  accurate  measurement 
procedures.  Secondly,  a  microchannel  heat  exchanger,  designed  and  built  specifically  for  this 
project,  was  tested.  In  the  context  of  this  study,  a  microchannel  heat  exchanger  is  defined  as  a 
heat  exchanger  with  small  primary  fluid  pore  spacing  and  thin  (normal  to  the  flow  direction) 
secondary  fluid  tubing.  A  detailed  description  of  both  heat  exchangers  follows. 

2.1  Straight  Plate  Finned-Tube  Heat  Exchanger 

A  straight  plate  fin-tube  heat  exchanger  was  built  by  a  commercial  manufacturer  and  is 
depicted  in  figures  land  2  along  with  relevant  dimensions.  This  heat  exchanger  was  part  of  a  set 
of  two  heat  exchangers  designed  and  built  for  a  thermoacoustic  prototype.  Due  to  manufacturing 
constraints,  it  was  necessary  to  make  several  compromises  in  the  design.  In  particular,  the 
stream  wise  depth  was  several  times  larger  than  the  acoustic  particle  displacement  associated 
with  the  thermoacoustic  prototype.  It  has  been  reported  that  depths  on  the  order  of  the  particle 
displacement  length  are  preferable  (Swift,  1999).  Furthermore,  the  width  and  spacing  of  the  fins 
were  not  optimal  according  to  guidelines  established  in  previous  design  optimization  studies 
(Minner,  1996). 
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Figure  2:  Picture  of  Installed  Fin-Tube  Heat  Exchanger 
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2.2  MicroChannel  Heat  Exchanger 

Figures  3  and  4  depict  the  microchannel  heat  exchanger.  Several  advantageous  features 
distinguish  this  design  from  that  of  the  conventional  fin-tube  design  shown  in  figure  1.  Firstly, 
the  use  of  small  rectangular  flow  channels  for  the  secondary  fluid  presents  a  smaller  obstruction 
to  the  flow  of  the  primary  thermoacoustic  working  fluid  for  a  given  tube  surface  area.  Secondly, 
the  use  of  the  microchannel  and  two-pass  secondary  flow  geometry  leads  to  higher  secondary 
fluid  heat  transfer  coefficients.  Thirdly,  the  microchannel  design  incorporates  outer  fins  having  a 
smaller  depth,  tighter  spacing,  and  larger  overall  surface  area. 


Figure  3:  Drawing  of  MicroChannel  Heat  Exchanger  (dimensions  in  inches) 
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Figure  4:  Picture  of  Installed  MicroChannel  Heat  Exchanger 


3.0  Experimental  Apparatus 

3.1  Wind  Tunnel  for  steady  flow  measurements 

A  small-scale  wind  tunnel  was  designed  and  built  to  perform  measurements  necessary  for 
estimating  air-side  heat  transfer  coefficients.  The  heat  exchangers  were  operated  at  steady  state 
in  a  steady  flow  with  water  flowing  through  the  tubes.  The  wind  tunnel,  shown  in  figures  5  and 
6,  consisted  of  a  centrifugal  fan  powered  by  a  Baldor  Electric  Motor  (Model  M31 15)  and  a  15 
cm  (6  in)  diameter  PVC  pipe.  A  flow  straightener  was  used  to  minimize  inflow  swirl  and 
transverse  flow  motion.  An  inlet  bell  mouth  and  a  plenum  settling  chamber  downstream  of  the 
heat  exchangers  were  utilized  to  improve  the  flow  velocity  profile  uniformity  across  the  duct. 

The  heat  exchanger  was  mounted  between  two  sections  of  PVC  pipe.  A  pitot  tube,  in 
conjunction  with  a  probe  traversing  mechanism  and  a  Dwyer  inclined  manometer  (Model 
number  115-AV),  was  used  to  measure  the  airflow  velocity  distribution  across  the  pipe  at  one 
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stream  wise  location.  The  velocity  distribution  was  used  to  determine  the  air  mass  flow  rate. 
Hot  water  flow  rates  through  the  heat  exchanger  were  measured  using  an  axial  paddle  wheel 
turbine  type  flow  meter  (JLC  International  IR-Opflow  Type  4)  at  the  water  outlet. 


Settling 

Chamber 


Water  in 


Air  Flow  Pitot  Tube 
Measurement 


Heat  Exchanger 


Flow 

Straightener 


Air  Outlet 
Temperature 
Measurement 


Water  out 


Air  Inlet 
Temperature 
Measurement 


Figure  5:  Schematic  of  the  Wind  Tunnel  used  for  the  Steady  Flow  Measurements. 


Figure  6:  Picture  of  the  Wind  Tunnel  used  for  the  Steady  Flow  Measurements. 

A  differential  temperature  transducer  (Delta-T  Company,  Model  number  75X)  was  used 
to  measure  the  water  temperature  difference.  At  the  inlet  and  the  outlet  of  the  wind  tunnel, 
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arrays  of  type  T  thermocouples  (N=9)  were  installed,  at  the  locations  shown  in  figure  5  to 
measure  air  temperatures.  The  rate  of  heat  transfer  was  calculated  from  the  measured  flow  rates 
and  temperatures  for  both  the  air-side  and  the  water-side  of  the  exchanger,  as  described  in 
section  5. 

3.2  Thermoacoustic  System  for  Unsteady  Flow  Measurements 

Unsteady  heat  exchanger  performance  measurements  were  made  using  a  functional 
thermoacoustic  cooler.  The  thermoacoustic  cooler  was  designed  and  built  at  Purdue  University 
(Minner  et  al.,  2000).  The  schematic  of  the  thermoacoustic  cooler  is  shown  in  figure  7  and  the 
picture  of  the  experimental  setup  is  shown  in  figure  8.  It  is  driven  by  a  300  W  moving  magnet 
linear  actuator  mounted  on  metal  "leaf"  springs  to  provide  suspension  stiffness.  Type  T 
thermocouples  were  used  to  perform  detailed  temperature  measurements  at  several  locations  on 
the  heat  exchangers.  The  driver  was  instrumented  with  an  accelerometer  on  the  driver  piston.  A 
pressure  sensor  was  installed  in  a  port  near  the  piston.  These  dynamic  pressure  and  velocity 
measurements  at  the  piston  face  were  used  as  inputs  to  DELTAE  (Ward,  1993)  to  estimate  the 
root-mean-square  (rms)  velocity  within  each  heat  exchanger.  This  velocity  was  used  to 
determine  acoustic  Reynolds  numbers  in  an  attempt  to  correlate  heat  transfer  data  for  periodic 
flows  and  to  compare  them  with  steady  flow  heat  transfer  coefficients,  as  described  in  greater 
detail  in  the  next  section. 

The  device  was  operated  with  two  heat  exchangers  and  a  stack  producing  a  temperature 
difference  between  the  circulating  water  stream  and  the  thermoacoustic  working  fluid  at  each 
heat  exchanger.  The  stack  was  composed  of  a  76  pm  thick  polyester  film  and  254  pm  thick 
nylon  wire  eonstructed  by  using  the  wire  as  a  spacer  and  rolling  the  film  into  a  cylinder.  The 
heat  transfer  rate  was  estimated  from  an  energy  balance  on  the  water  stream  using  measurements 
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obtained  with  the  same  instrumentation  described  in  the  previous  section.  The  temperatures  at 
the  end  of  the  stack  adjacent  to  each  heat  exchanger  were  measured  with  thermocouples  located 
near  the  centerline. 


Hot  Heat 

Moving  Coil  Exchanger 
Transducer 


Cold  Heat 
Exchanger 


Piston 


t 


Closed  End 


Stack 


Figure  7:  Schematic  of  Thermoacoustic  Cooler. 


Figure  8:  Picture  of  the  Thermoacoustic  Cooler. 
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4.0  Experimental  Procedures 


4.1  Steady  Flow  Experimental  Procedures 

The  steady  flow  tests  were  designed  to  conform  to  the  ARI  Standard  410  for  Forced 
Circulation  Air-Cooling  and  Air-Heating  Coils  (www.ari.org).  This  standard  establishes  the 
proper  flow  conditions  of  the  primary  and  secondary  loops,  and  outlines  procedures  for  varying 
the  flow  conditions. 

The  data  were  gathered  using  a  Hewlett  Packard  75000  Series  B  data  acquisition  system 
with  inputs  from  the  two  thermocouple  arrays,  the  differential  temperature  transducer,  the  water 
flow  meter,  and  eight  thermocouples  distributed  over  the  heat  exchanger.  Hewlett  Packard’s 
visual  design  environment  HPVEE  was  used  to  control  the  data  acquisition.  Each  data  point  was 
marked  with  a  time  reference  and  was  stored  in  a  text  file  associated  with  a  particular  data  set. 
The  data  was  then  post-processed  using  an  Excel  spreadsheet. 

4.2  Oscillating  Flow  Experimental  Procedures 

Data  acquisition  procedures  in  the  thermoacoustic  environment  were  similar  to  those 
followed  in  the  steady  flow  environment.  The  heat  transfer  rate  was  calculated  from  an  energy 
balance  and  water  side  temperature  data.  It  was  necessary  to  use  measured  dynamic  pressure  and 
piston  velocity  data  to  determine  the  particle  velocity  in  the  heat  exchangers  as  discussed  later. 
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5.0  Estimation  of  Heat  Transfer  Coefficients 


For  heat  transfer  in  steady  flow,  a  local  heat  transfer  coefficient  is  defined  as  the  ratio  of 
the  heat  transfer  per  unit  area  to  the  temperature  difference  between  the  surface  and  “bulk”  fluid 
adjacent  to  the  surface.  However,  in  practice,  a  global  heat  transfer  coefficient  is  utilized  that  is 
derived  from  a  lumped  analysis  of  the  heat  exchanger  (UA-LMTD  or  effectiveness-NTU 
methods).  Generally,  a  global  heat  transfer  coefficient  is  estimated  from  an  overall  heat 
exchanger  conductance  (UA)  determined  from  measurements,  using  existing  correlations  for 
secondary  fluid  heat  transfer  coefficient  and  overall  fin  efficiency.  The  UA  of  the  heat 
exchanger  is  estimated  from  the  measurements  as  the  ratio  of  heat  transfer  rate  to  a  log-mean 
temperature  difference  between  the  primary  and  secondary  fluids.  The  overall  heat  transfer  rate 
is  estimated  from  primary  or  secondary  fluid  temperature  differences  and  flow  rate 
measurements.  The  heat  transfer  coefficient  data  are  often  non-dimensionalized  using  the  j- 
Colbum  factor.  This  non-dimensional  presentation  allows  extrapolation  to  different  gases  and 
operating  conditions.  The  j-Colbum  factor  data  are  usually  correlated  in  terms  of  the  Reynolds 
number. 

For  the  purpose  of  modeling  heat  exchangers  within  a  thermoacoustic  device,  it  is  more 
convenient  to  define  a  local  heat  transfer  coefficient  as  the  ratio  of  the  heat  transfer  per  unit 
surface  area  to  the  temperature  difference  between  the  heat  exchanger  surface  and  the  time- 
average  thermoacoustic  fluid  temperature  at  the  end  of  the  stack  adjacent  to  the  heat  exchanger. 
Previous  researchers  have  estimated  heat  transfer  coefficients  for  thermoacoustic  heat 
exchangers  by  assuming  conduction  across  a  laminar  boundary  layer  (Minner,  1996).  However, 
results  from  these  calculations  have  never  been  compared  with  measured  values. 


21 


In  the  current  study,  global  thermoacoustic  fluid  heat  transfer  coefficients  were  estimated 
from  overall  heat  exchanger  conductances  (UA)  determined  from  measurements,  using  existing 
correlations  for  secondary  fluid  heat  transfer  coefficient  and  overall  fm  efficiency.  This  global 
heat  transfer  coefficient  is  consistent  with  a  lumped  analysis  of  the  heat  exchanger  that  includes 
the  effects  of  temperature  variations  within  the  fins  and  secondary  fluid.  The  UA  of  the  heat 
exchanger  was  estimated  from  the  measurements  as  the  ratio  of  heat  transfer  rate  to  a  log-mean 
temperature  difference  between  the  secondary  fluid  and  the  time-averaged  thermoacoustic  fluid 
temperature  near  the  stack.  The  overall  heat  transfer  rate  was  estimated  from  secondary  fluid 
temperature  and  flow  rate  measurements.  The  heat  transfer  coefficient  data  were  non- 
dimensionalized  using  the  j-Colbum  factor.  In  this  case,  the  j-Colbum  factor  data  were 
correlated  in  terms  of  an  acoustic  Reynolds,  defined  using  a  root-mean-square  (rms)  velocity 
within  the  heat  exchanger.  The  rms  velocities  at  the  heat  exchangers  were  estimated  from  linear- 
acoustic  theory  using  dynamic  pressure  and  velocity  measurements  at  the  driver  face  of  the 
thermoacoustic  device  and  DELTAE  (Ward,  1993). 

The  following  subsections  provide  details  of  the  calculations  used  to  estimate  heat 
transfer  coefficients  from  both  physical  descriptions  and  measured  performance  data. 

5.1  Steady  Flow  Calculations 

5.1.1  Theoretical  Calculations 

Many  general  correlations  exist  for  estimating  global  heat  transfer  coefficients  for  steady 
flow  over  fin-tube  geometries  given  physical  dimensions  and  properties.  For  heat  exchangers 
with  smooth  fins,  a  correlation  developed  by  Gray  and  Webb  [1986]  expresses  heat  transfer  data 
in  terms  of  a  j-Colbura  factor  according  to  the  following  correlation. 


22 


Nu 


RePr'^^ 


.991 X 


2.24Re"°’2 


/jy.03A 

.4. 


1.82  A 


.14  Re 


-.328 


-.502 


v^oy 


.0312A 


(1) 


where  the  tube  spacing  normal  to  the  flow  is  given  by  St,  the  tube  spacing  in  the  air  flow 
direction  is  given  by  St  (This  factor  was  ignored  in  this  study  because  there  is  one  row  in  the  heat 
exchanger),  the  spacing  between  adjacent  fins  is  s,  D„  represents  the  tube  diameter  for  the 
secondary  flow,  and  Nu,  Re,  and  Pr  refer  to  standard  dimensionless  groups  defined  in  the 
nomenclature. 


5.1.2  Calculations  from  Data 

The  overall  conductance  of  a  heat  exchanger  operating  in  steady  flow  and  at  steady  state 
can  be  estimated  from  experimental  data  using  a  log-mean  temperature  difference  model.  A 
cross-flow  heat  exchanger  is  considered  in  this  study  with  water  flowing  through  the  tubes  and 
air  flowing  over  the  finned  tubes.  For  the  operating  conditions  considered  in  this  study,  the 
capacitance  rate  (product  of  mass  flow  rate  and  specific  heat)  of  the  water  is  much  greater  than 
that  of  the  air.  In  this  case,  the  overall  heat  exchanger  conductance  is  estimated  according  to 


UA  = 


Q 

LMTD 


(2) 


where  Q  is  the  heat  transfer  rate  determined  from  measurements  and  LMTD  is  the  log-mean 
temperature  difference  defined  as 


LMTD, 


In 


T  -T 


(3) 


and  where  s  denotes  steady  flow  and  7),„  Td,  Tho,  and  Tco  are  the  inlet  (i)  and  outlet  (o) 
temperatures  of  the  hot  Qi)  and  cold  (c)  streams. 
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The  heat  transfer  rate  is  determined  from  measurements  using  an  energy  balance  on  the 
hot  or  cold  stream  as 

Q  =  ^hCph(TH.i-Th.„)  (4) 

or 

Q  =  th,cJT^„-T^i)  (5) 

where  m  is  mass  flow  rate,  Cp  is  constant  pressure  specific  heat,  and  the  subscripts  h  and  c 
denote  hot  or  cold  stream  conditions,  respectively. 

The  heat  exchanger  conductance  can  be  estimated  from  measurements  using  equations  2, 
3,  and  4  (or  5).  The  heat  transfer  conductance  can  then  be  related  to  the  hot  and  cold  stream  heat 
transfer  characteristics  according  to 


UA  =  - 


■ + +  - 


hA 


(6) 


where  t],  h,  and  A  are  overall  fin  efficiency,  global  heat  transfer  coefficient,  and  surface  area,  Rtw 
is  the  resistance  of  the  tube  wall,  and  the  subscripts  o  and  i  denote  outer  and  inner  surfaces  of  the 
fin-tube  heat  exchanger. 

Generally,  the  outer  global  heat  transfer  coefficient  {hg)  is  estimated  from  data  by  solving 
equation  6  using  a  UA  determined  from  measurements,  known  surface  areas,  and  theoretical 
estimates  for  Rtw,  hi,  and  r]o. 
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The  inside  heat  transfer  coefficient,  hi,  for  the  fin  tube  heat  exchangers  can  be 
conveniently  found  for  a  range  of  Reynolds  number  (3000<  ReD<5  x  10®)  using  the  Gnielinski 
correlation  (Incropera  and  Dewitt,  1996) 


(f^/8)(  Rep-1000)Pr 

l  +  12.7(/c/8)'"VPr^'^-l) 


(7) 


where  Nud  is  the  Nusselt  number  and  is  related  to  the  heat  transfer  coefficient  according  to 

(8) 

k 

where  hi  is  the  inside  convection  heat  transfer  coefficient,  k  is  the  thermal  conductivity  and  D  is 
the  inner  diameter  of  the  finned  tube.  Therefore  hi  can  be  calculated  from  Egs.  7  and  8. 

Also,  the  friction  factor, /g,  can  be  found  by 

/c=f.791nRe„-1.64)-^  (9) 

and  the  Reynolds  number  based  up  tube  diameter  is  given  by 

(JO, 

where  u  is  the  fluid  velocity,  D  is  the  inner  diameter  of  the  finned  tube,  p  is  the  fluid  density,  and 
//  is  the  dynamic  viscosity.  Due  to  the  different  tube  geometry  of  the  microchannel  design, 
Modine  Manufacturing  Corporation  recommends  using  the  Sieder-Tate  correlation. 


Nud  =1.861^^^^’^ 


where  L  is  the  tube  length,  //  is  the  dynamic  fluid  viscosity,  is  dynamic  viscosity  of  the  fluid  at 
the  tube  surface,  and  Pr  is  the  Prandtl  number  and  D/,  is  the  hydraulic  diameter  of  the  non 
circular  tube  which  is  defined  as 
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(12) 


4A, 


where  A„  is  the  cross-sectional  area  and  P  is  the  wetted  perimeter  of  the  finned  tube. 

The  tube  wall  resistance  and  external  fin  efficiency  are  basic  quantities  found  in  most 
heat  transfer  textbooks  (Incropera  and  Dewitt,  1996)  and  are  defined  by: 


27tL,k, 


(13) 


and 


(14) 


where 


tanh(mL) 

mL 


(15) 


and 


m  = 


(16) 


In  equations  13-16,  Do  and  D,  are  the  external  and  internal  diameters,  L,is  the  total  tube  length,  ife, 
is  the  thermal  conductivity  of  the  tube.  A/ is  the  total  fin  surface  area,  A  is  the  total  heat 
exchanger  surface  area,  L  is  the  length  of  a  single  fin,  P  is  the  fin  perimeter,  ho  is  the  external 
heat  transfer  coefficient,  k  is  the  thermal  conductivity  of  the  fin  material,  and  A/c  is  the  fin  cross- 
sectional  area. 

In  order  to  non-dimensionalize  the  data  for  various  gases  and  to  facilitate  comparison 
with  theoretical  calculations,  the  experimental  outer  heat  transfer  coefficient  are  converted  into 
the  j -Colburn  factor  define  as: 
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(17) 
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where  Nud  and  Re^,  are  Nusselt  number  and  Reynolds  number,  respectively,  for  the  air  side 
flow.  The  hydraulic  diameter,  ,in  this  case  is  conveniently  calculated  by 


D 


h 


4LA, 

A 


(18) 


where  L  is  the  fin  length  (length  of  flow  path),  Ac  is  the  flow  area  and  A  is  the  total  heat 
exchanger  surface  area.(Gray  and  Webb,  1986).  The  detailed  calculation  for  D/,  is  given  in 
appendix  IX 

5.2  Description  of  Oscillating  Flow  Calculations 
5.2.1  Theoretical  Calculations 

A  general  definition  for  local  heat  transfer  coefficient  for  a  surface  in  contact  with  a  fluid 
in  internal  flow  is: 


-kdT/dy\  0 


(19) 


where  k  is  the  thermal  conductivity  of  the  fluid,  Ts  is  the  surface  temperature  and  r„,  is  the  mean 
temperature  of  the  fluid  over  the  surface  of  the  heat  exchanger.  In  order  to  provide  a  theoretical 
method  for  calculating  heat  transfer  coefficients  in  an  oscillating  flow,  several  investigators  have 
proposed  the  use  of  equation  19  with  the  mean  fluid  temperature  taken  to  be  the  time-averaged 

fluid  temperature  at  the  stack  end  adjacent  to  the  heat  exchanger,  3T/3y|  y^o  approximated  by 

the  temperature  difference  between  the  surface  and  the  time-averaged  stack  end  fluid 
temperature  divided  by  the  thermal  boundary  layer  thickness,  (Minner,  1996).  The 
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boundary  layer  thickness  is  assumed  to  be  the  minimum  of  the  half  plate  spacing  of  the  heat 
exchanger  (yo)  and  the  thermal  penetration  depth  (<5{)  where 


With  these  approximations,  equation  17  reduces  to 


rc\m(y„,6, ) 

For  an  ideal  thermoacoustic  heat  exchanger,  only  the  surface  area  within  one  particle 
displacement  of  the  stack  can  participate  in  the  heat  transfer.  Therefore,  a  global  heat  transfer 
coefficient  that  accounts  for  the  “dead”  surface  area  of  the  heat  exchanger  is  estimated  as 


minry„,^J 

where  Xp  is  the  distance  associated  with  a  particle  displacement  at  the  heat  exchanger  and  is 
the  width  of  the  heat  exchanger  in  the  direction  of  thermoacoustic  flow  (Minner,  1996). 

It  has  also  been  proposed  that  heat  transfer  coefficients  for  oscillating  flow  conditions 
could  be  estimated  using  steady-flow  correlations  with  the  Reynolds  number  evaluated  using  the 
acoustic  velocity  (rms  of  the  oscillating  velocity)  (Mozurkewich  (1995).  Both  the  boundary 
layer  approximation  and  use  of  steady-flow  correlations  are  investigated  as  approaches  for 
predicting  heat  transfer  coefficients  in  oscillating  flow. 


28 


5.2.2  Calculations  from  Data 


For  cases  where  the  heat  exchanger  is  placed  in  an  acoustically  oscillating  flow  (within  a 
thermoacoustic  system),  equations  2,  3, 4  (or  5)  and  6  are  also  used  to  estimate  the  outer 
(thermoacoustic  fluid  to  surface)  heat  transfer  coefficient  from  measurements  and  the  heat 
transfer  coefficient  is  transformed  into  the  j-Colbum  factor  using  equation  17.  However,  the 
inlet  and  outlet  primary  fluid  temperatures  are  not  defined  for  the  oscillating  flow  stream  and 
little  temperature  change  can  be  attributed  to  fluid  particles  oscillating  either  over  the  heat 
exchanger  or  between  the  heat  exchanger  and  the  stack.  In  this  case,  the  log  mean  temperature 
difference  was  calculated  by 


LMTD,^^ 


(T  -T  ■) 

'  w,o  w,i  / 


T  -T 

Lni-^ - 

\Ts-T.,o\ 


(24) 


where  tac  stands  for  thermoacoustic,  Tw.o  and  Tw,i  are  the  water  outlet  and  inlet  temperatures,  and 
Ts  is  the  stack  end  temperature.  Finally,  the  steady-state  energy  balance  for  the  water  stream  is 
used  to  determine  the  heat  transfer  rate  from  measurements. 

In  this  study,  root  mean  square  (rms)  acoustic  particle  velocity  was  used  to  correlate  j- 
Colbum  factors  for  oscillating  flow  conditions.  Unfortunately,  it  was  not  possible  to  measure  the 
acoustic  particle  velocity  in  the  heat  exchangers  directly.  This  problem  was  solved  by  using  a 
thermoacoustic  modeling  software  developed  at  Los  Alamos  National  Laboratory,  DeltaE  (Ward 
1993),  to  compute  the  local  root-mean  square  velocity.  This  estimate  of  the  particle  velocity  was 
then  used  in  computing  an  acoustic  Reynolds  number.  The  acoustic  Reynolds  number  was  used 
in  correlating  the  j-Colbum  factors. 
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6.0  Results 


6.1  Steady  Flow  Results 

The  data  for  the  steady  flow  results  are  presented  in  two  forms.  Figure  9  is  a  graph  of  the 
j-Colbum  factor  versus  the  Reynolds  number.  Figure  10  shows  the  heat  transfer  coefficient 
versus  the  Reynolds  number.  In  each  case,  the  heat  transfer  coefficients  were  calculated  using 
heat  transfer  results  determined  from  air-side  and  water-side  measurements.  For  each  airflow 
rate,  several  different  water  flowrates  were  used  in  order  to  verify  consistency.  The  differences 
between  the  air-side  and  water-side  results  are  an  indication  of  experimental  errors.  Generally, 
the  air-side  and  water-side  results  are  within  about  20%. 

The  heat  transfer  coefficient  increases  with  Reynolds  number,  as  expected.  The  fin-tube 
heat  exchanger  results  were  compared  with  the  Gray-Webb  correlation,  shown  as  a  solid  line  in 
figure  9.  The  test  results  compare  very  well  with  the  correlation.  No  similar  correlation  could  be 
found  in  the  literature  for  the  micro-channel  heat  exchanger.  The  heat  transfer  coefficients  for 
the  microchannel  heat  exchanger  are  significantly  greater  than  those  of  the  fin-tube  heat 
exchanger.  Note  the  different  hydraulic  diameter  for  the  two  exchangers. 
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Figure  9.  Steady  Flow  j-Colbum  v^.  Reynolds  Number 


Reynolds  Number 


Figure  10.  Steady  Flow  Heat  Transfer  Coefficient  (W/m^'K)  v^.  Reynolds  Number 
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6.2  Performance  of  the  Heat  Exchangers  in  Oscillating  Flow 

The  results  for  the  fin-tube  heat  exchanger  only  are  shown  in  figures  1 1  and  12,  where 
the  j-factor  and  the  heat  transfer  coefficient  h  are  plotted  versus  acoustic  Reynolds  number, 
respectively.  Unfortunately,  heat  gains  originating  from  the  driver  biased  much  of  the  early  data 
and  only  a  very  small  subset  of  the  original  data  was  useable.  In  order  to  determine  which  data 
were  unusable,  a  first  law  balance  was  used.  The  sum  of  the  measured  power  input  and  the 
pumped  heat  were  compared  to  the  measured  heat  rejected  on  the  hot  side  heat  exchanger.  Data 
points  with  an  energy  balance  disparity  greater  than  15%  were  rejected.  As  can  be  seen  from  the 
fin-tube  results,  only  a  limited  number  of  data  points  passed  the  test.  This  caused  poor  resolution 
in  the  data.  Nevertheless,  the  data  clearly  shows  a  marked  increase  in  performance  for  the 
microchannel  heat  exchanger.  Also,  it  is  obvious  that  the  boundary  layer  approximation 
provides  estimates  for  heat  transfer  coefficients  that  are  extremely  conservative. 
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Figure  II.  j-Colbum  V5.  Reynolds  Number  for  Periodic  Flow;  Fin-Tube  Heat  Exchanger. 
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Figure  12.  Heat  Transfer  Coefficient  (W/m^'K)  vs.  Reynolds  Number  for  Periodic  Flow;  Fin- 

Tube  Heat  Exchanger 


Figures  13  and  14  show  the  results  for  the  micro-channel  heat  exchangers.  In  each  case, 
results  for  both  the  hot  and  cold  heat  exchangers  in  the  thermoacoustic  system  are  presented.  In 
addition  to  the  experimental  results,  the  figures  show  the  boundary  layer  conduction 
approximation  that  is  often  used  in  thermoacoustics  to  predict  heat  transfer.  For  the 
microchannel  heat  exchanger,  several  tests  were  performed  where  heat  leak  from  the  driver  was 
kept  relatively  small.  In  order  to  reduce  the  effects  of  heat  gains  from  the  driver,  the  driver  was 
only  operated  short  periods  of  time  for  each  data  point.  Using  this  method,  it  was  possible  to 
reduce,  but  not  completely  eliminate,  the  energy  imbalances.  As  can  be  seen  from  figures  13  and 
14,  the  data  appears  to  be  consistent,  and  several  interesting  trends  appear.  Firstly,  the  measured 
heat  exchanger  performance  at  low  Reynolds  numbers  is  much  greater  than  that  predicted  using 
the  boundary  layer  calculation  commonly  used  in  thermoacoustics.  As  mentioned  earlier,  these 
models  assume  that  the  heat  transfer  coefficient  is  a  result  of  conduction  across  the  thermal 
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boundary  layer.  It  appears  that  this  approach  works  well  at  high  Reynolds  numbers,  but 
significantly  underpredicts  the  j-Colbum  factor  at  low  Reynolds  numbers.  This  could  have  a 
significant  impact  on  the  prediction  and  the  design  optimization  of  the  performance  of 
thermoacoustic  systems. 

A  second  observation  is  that  the  steady  flow  regression  model  seems  to  provide  a 
reasonably  good  approximation  for  the  thermoacoustic  j-Colbum  factor  when  evaluated  with 
acoustic  Reynolds  numbers.  It  is  obvious  from  the  figure  that  the  steady  flow  regression 
provides  a  better  estimate  than  the  boundary  layer  approximation.  This  can  be  also  observed  in 
figure  15 which  compares  the  regression  models  for  the  thermoacoustic  and  steady  flow  data  of 
the  microchannel  heat  exchanger. 

Another  important  trend  is  a  difference  in  the  functional  dependence  on  Reynolds  number 
between  the  hot  and  cold  heat  exchangers.  It  is  believed  that  this  difference  is  due  an  inaccurate 
characterization  of  the  stack-end  temperatures.  In  order  to  investigate  this  further,  a  radial  array 
of  thermocouples  was  installed  at  the  stack  ends.  This  array  revealed  that  there  is  an  obvious 
dependence  of  stack  temperature  upon  radial  position.  It  has  been  postulated  that  this  is  due  to 
the  influence  of  the  container’s  aluminum  walls  and  it  obviously  implies  that  using  one 
temperature  probe  to  calculate  the  log-mean  temperature  difference  is  misleading  since  it  does 
not  properly  take  into  account  the  temperature  distribution.  Using  data  from  the  thermocouple 
array,  a  spatially  averaged  temperature  distribution  was  calculated  and  used  to  calculate  a  new  j- 
Colbum  factor.  The  results  appear  in  figure  16  for  a  limited  number  of  data  points.  It  can  be 
seen  that  this  new  approach  yields  closer  agreement  between  the  hot  and  cold  heat  exchanger 
heat  transfer  coefficient  relationships.  However,  in  this  case,  the  steady-flow  correlation  tends  to 
over  predict  the  heat  transfer  coefficients. 
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Colburn  Factor 


Finally,  the  data  show  that  the  microchannel  heat  exchanger  performs  significantly  better 
than  the  fin  tube  heat  exchanger  when  operating  in  an  oscillating  flow  environment.  These 
results  are  shown  in  Figure  17  and  are  consistent  with  the  steady-flow  results. 


Reynolds  Number 

Figure  13.  j-Colbum  v^.  Reynolds  Number  for  Periodic  Flow;  MicroChannel 
;  Dotted  line  denotes  the  extrapolation  of  the  experimental  data 
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i-Colbum  Factor  ‘’l  Heat  Transfer  Coefficient 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000  10000 

Reynolds  Number 


Figure  15.  j-Colbum  vs.  Acoustic  Reynolds  Number  Regression;  MicroChannel  Heat  Exchanger 
; Dotted  line  denotes  the  extrapolation  of  the  experimental  data 


36 


j-Colbum  Factor  (#5"  j-Colbum  Factor 


ire  16.  j-Colbum  vs.  Acoustic  Reynolds  Number;  MicroChannel  Heat  Exchanger,  Computed 

from  Radial  Thermocouple  Array  Data 
; Dotted  line  denotes  the  extrapolation  of  the  experimental  data 


0  1000  2000  3000  4000  5000  6000  7000  8000  9000  10000 

Reynolds  Number 


Figure  17.  j-Colbum  V5.  Acoustic  Reynolds  Number:  Comparison  of  Heat  Exchangers 
; Dotted  line  denotes  the  extrapolation  of  the  experimental  data 
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7.0  Conclusions  and  Future  Work 


The  performance  of  two  different  heat  exchangers  was  measured  in  steady  and  pulsating 
flows.  The  results  suggest  that  the  thermoacoustic  heat  transfer  could  be  predicted  using  steady 
flow  correlations  with  acoustic  Reynolds  numbers.  This  could  simplify  the  heat  transfer 
calculations  and  would  increase  the  accuracy  of  the  current  models.  In  order  to  confirm  this,  a 
new  battery  of  tests  would  need  to  be  designed  and  implemented.  It  would  be  necessary  to  run 
the  tests  on  several  different  heat  exchanger  types  and  compare  them  to  the  steady  flow 
correlations.  It  should  be  noted  that  there  is  an  uncertainty  in  the  value  of  the  particle  velocity 
because  it  is  being  estimated  from  linear  acoustic  models  (DeltaE  calculation).  This  value  could 
also  be  affected  by  possible  streaming  (i.e.  acoustically  driven  steady  flows)  effects  within  the 
thermoacoustic  device.  These  possible  errors  could  be  assessed  and  corrected  by  implementing  a 
hot-wire  anemometer  into  the  system  at  specific  points  to  measure  the  total  (AC  and  DC) 
velocity. 

The  accuracy  of  the  test  data  for  the  pulsating  flow  cases  is  strongly  dependent  on  the 
magnitude  of  the  heat  flux.  In  its  current  configuration,  the  thermoacoustic  system  does  not 
allow  very  high  heat  fluxes  to  be  achieved,  the  the  system  is  affected  by  heat  production  from  the 
driver,  and  thermal  conduction  though  the  shell.  Temperatures  gradients  are  small.  It  is  difficult 
to  accurately  measure  small  temperatures,  especially  when  considering  that  the  thermoacoustic 
system  is  closed  and  pressurized.  The  best  solution  for  eliminating  this  source  of  error  is  to 
continue  to  work  to  increase  the  heat  transfer  values  and  install  a  driver  cooling  system. 
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APPENDICES 


Appendix  I:  Table  of  Thermoacoustic  Raw  Data  for  the  MicroChannel  Heat  Exchangers 


Measurement  [units 


Wac 


Qc 


Qh 


DTh 


DTc 


Vrateh 


Vratec 


stack  DT 


0 


Qc+Wac-Qh 


hi 


wtr-stk 


wtr-stk 


iston  vei 


r  Chase 
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Mtbh 


Mtbc 


inbai/(Qc+Wac)* 
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Mfn-h 
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GBFSC-c 
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CO 
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Whi 


Wei 


IdegC 


deoC 


DT  -h  dea  C 


DT  -c  dea  C 


200 


0.55 


172.10 


65.77 


91.17 


144.91 


0.73 


0.56 


0.05 


0.04 


5.26 


53704.05 


12.03 


22.91 


22.81 


-2.91 


2.26 


0.96 


125.03 


23.28 
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25.73 


19.77 


7.67 


26.16 


20.57 
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22.25 
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22.81 


3 


200 


0.55 


172.1 


36.87 


62.00 


95.14 


0.481  1.151 


Data  Point 


6  7 


300  300 


0.55  0.55 


173.25  173.25 


IjcMvMIRKTlKHSISIBjgBtli 


0.38 


0.05 


0.04 


3.49  7.47 


0.98  1.12 


0.54  0.61 


0.05  0.05 


0.04  0.04 


6.16  6.80 


31417.2 


3.73 


22.93 


22.81 


-1.84! 


1.54 


0.72 


124.16' 


23.17 


-2.46  0.89 


22.97 


22.84  22.83 


-3.61  -4.03 


2.43  2.64 


1.07  1.15 


UliWI 


19.90 


27.28 


20.39 


26.57 


20.41 


23.97  24.14 
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24.28  24.49 


22.97 


8 


300 


0.3447 


143.90 


83.71 


84.94 


179.59 


0.90 


0.50 


0.05 


0.04 


5.71 


118582.5 


-10.94 


22.88 


22.75 


-3.30 


2.28 


0.97 


-119.03 


23.33 


23.00 
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6.49 


26.63 


20.69 


19.81 


26.53 


20.38 


27.04 
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26.18 


20.48 


23.85 
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22.75 


9 
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119.14 
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0.60 
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6.95 
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-5.32 
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2.71 


1.11 
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23.45 


23.06 


26.83 


19.62 


2.40! 


27.42 


20.40 


19.34 


27.25 


20.05 


28.19 


20.22 


27.00 


20.05 


24.10 


22.181 


20.33 


24.52 


22.88 


22.76 


Measurement 


Qc 


Qh 


DTh 


DTc 


Vrateh 


Vratec 


stack  DT 


Hz 


W 


W 


W 


deaC 


degC 


L/sec 


L/sec 


degC 


kPa 
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C 


C 


13 


300 


0.216 


133.75 


80.99 


75.24 


167.07 


0.84 


0.42 


0.05 


0.04 


5.10 


121013.5 


-10.84 


22.94 


22.82 


14 


300 


0.216 


133.85 


116.82 


87.31 


200.24 


1.01 


0.48 


0.05 


0.04 


6.10 


146558.5 


3.89 


22.94 


22.82! 


16 


250 


0.216 


133.70 


88.47 


85.32 


186.78 


0.94 


0.47 


0.05 


0.04 


5.74 


107562.5 


-12.99 


23  01 


22.90 


Data  Point 


18l  19 


200 


0.55  0.55 


168.50 


110.43 


116.78 


246.52 


1.16  1.25 


0.64 


0.05  0.05 


0.04 


7.20 


79362.5 


-19.30 


23.01  23.00 


22.89 


20 


292 


0.466 


162.00 


21.58 


39.43 


49.71 


0.24 


0.26 


0.05 


0.04 


1.94 


22155.4 


11.30 


22.55 


22.45! 


22424.7 


18.84 


22.61 


22.51 


22 
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0.55 


170.25 


41.30 


61.72 


112.28 


0.54 


0.40 


0.05 


0.04 


3.26 


19393.951 


23 


250 


0.55 


170.25 


26.38 


43.18 


78.88 


0.38 


0.28 


0.05 


0.04 


2.39 


18770.1 


24 


313 


0.481 


162.30 


55.51 


82.51 


122.99 


0.56 


0.51 


0.05 


0.04 


4.09 


84432.5 


15.03 


22.60 


22.50 
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GBFDC-c 


GBFSQ-h 


hi 


Wei 


21.58  21.24 


23.83  24.61 


22.72  22.77 


21.65  21.38  21.77 


24.07  24.72  23.92 


22.93  22.76  22.72 


22.83  22.65  22.63 


Measurement  lunits 


Mtbh 


Mtbc 


inbal/(Qc+Wac)' 


GBSI-h 


GBS2-h 


GBS3-h 


GBS1-C 


GBS2-C 


GBS3-C 


Hot  Stack  Av 


Cold  Stack  Av 


35 


200 


0.55 


169.70 


75.41 


91.72 


150.55 


0.90 


0.65 


0.04 


0.03 


5.49 


55968.05 


16.59 


24.02 


23.76 


-4.20 


2.09 


0.81 


140.77 


24.47 


24.09 


28.62 


28.53 


9.92 


27.77 


28.53 


28.62 


Data  Point 


36  37 
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0.55  0.55 
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75.96 
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173.90 


0.98  1.11 


0.66 


0.04 


0.03 


5.47 


55479.2 


-6.09 


24.00 


23.77 


-4.27 


2.10  2.44 


0.80 


142.22 


24.49  24.58 


24.10  24.17 


28.68 


28.56 


3.63 


27.851  28.11 


28.57 


28.68 


21.92 


21.671  21.33 


21.55  21.28 


28.27 


38 
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0.55 


169.70 


69.63 


90.23 


179.38 


1.01 


0.65 


0.04 


0.03 


5.45 


45165.7 


39 
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0.55 


169.90 


76.50 


97.93 


201 .30 


1.13 


0.71 


0.04 


0.03 


6.24 


44257.55 


-26.88 


24.06 


23.84 


-4.33 


2.43 


0.93 


139.52 


24.63 


24.19 


28.79 


28.69 


15.41 


28.02 


28.65 


28.79 


21.65 


21 .42 


21.32 


28.40 


21.41 


40 


200 


0.55 


169.50 


70.04 


92.27 


194.89 


1.10 


0.66 


0.04 


0.03 


6.19 


36761.65 


-32.59 


24.22 


23.98 


-3.90 


2.49 


Key  to  Notation 


D- temperature  difference 


o-out 

h-hot 

c-cold 

M-metal 

tb-tube 

f-fin 

G-gas 

B-between 

F-fins 

D-duct  side 


S -stack  side 

38-3/8”  off  of  the  heat  exchanger 
C-center  of  heat  exchanger 
Q-quarter  diameter  of  heat  exchanger 
Wac-acoustic  power 
Qc,  Qh  -  heat  transferred 
V-voIumetric  flow  rate 
pO-pressure  at  driver 

wtr-stk  -  temp  diff  between  water  and  stack 
GBSl-3  -  temp  of  gas  stack  for  area  ratio 
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Appendix  II:  Table  of  Thermoacoustic  Results  for  the  MicroChannel  Heat  Exchangers 


Data  Point 


1 


2 


3 


Hot  Exchanger  Results 


Acoustic  Information 


1727.93 


BHEBBESMltHBclIiag! 

BHIW^lticicmciriiKggl^fgcai 

BEHKEiHfeamsirfliSiaEai 

ECgg>SK51i^Mglli!Sfgi 

BSlEEElEEB^EgKgral 

lEEBBKSBEBKBUrfefcal 

BBligRilj^^E^TOKBcgj 

KEHin^ljcb^g^ia.ncwi 

liEHmajMKaasaKtS  11^1=1^1^11 

lhit:wWE5BB!SHro^lS5llcl.-ilicKii 

EElBEgggBiWTBtCTai 


and  Calculation 


h  air  side 


781 .61 


764.21 


879.83 


1091.46 


797.97 


1045.40 


1183.28 


1073.96 


1123.59 


1305.63 


1242.21 


1217.31 


1213.35 


1128.16 


1215.88 


1135.43 


1301.15 


1364.08 


air  side  |H20  vel  1 


0.0180 


Water  Flow  in  HX  Tubes  Calculation 


channel  iReD  Kriction  factor 


I^KS 


1493.17 

0.0269 

1.0251 

625.72 

0.0713 

1.0746 

1410.00 

0.1241 

1.0746 

1942.63 


1606  08 


1114.041 


mmi 

DBI 

mmi 

lEm 

IIEIESI 

IDSI 

\mmn 

ilTMl 

\mmi 

mmi 

DSI 

IKII 

irffiimi 

mmi 

mmi 

nm 

in^i 

mmi 

■Basil 

mil 

mmi 

mmi 


BBBll! 


1 .0746 


1.07461 1937.13 


1.1312 


1.1312 


1.1312 


1.1312 


1.1312 


1.1312 


1.1312 


1.1312 


0.8661 


0.8661 


0.8661 


0.8661 


0.91911  1731.79 


imm 

imxE 

SEm' 

i^m\ 

3EMI 

mm 

mm\ 

ima 

im 


msam 


3717.03 


3745.72 


lisiisia^ 

iiiT^BgnEg 


10.41 

EEeEE] 

10.37| 

13714.03 

0.0139 

0.9191 

1719.13 

0.0555 

0.0182 

0.9191 

■ECBBH 

0.0555 

ii^si 


Ml»g!jBg!SfB 

IMiKilE^E 

SMBgnianBg 

]EEM 

BtbkliEgEl 

B|bl^B:!c|:^4il 

|K^BKi;icnmi 

HEEIESES 

gMt^ESnCH 

M|^E!3Bnil 

Mii^BSHna 

Ml»>i!.-iESBEg 


Cold  Exchanger  Results 


_ Acoustic  Information  and  Calculation 


Water  Flow  in  HX  Tubes  Calculation 
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295.59  0.0041  0.9368  1705.57 


290.58  I  0.0039  0.9368  1708.12 


BBlIcHCTi 


friction  factor  Nusselt 


htubes 


3590.10 


3591.03 


3592.46 


3588.48 


3591.59 


3590.69 


3590.14 


3641.68 


3640.70 


3709.81 


3708.77 


44 


12 

EBl 

6974.17 

1  1.86 

34.11 

249.33 

0.0027 

0.9368 

1696.37 

0.0558 

9.92 

3713.57 

13 

ES! 

EBl 

7992.87 

E3!! 

253.50 

EESB 

0.9368 

1698.59 

0.0557 

9.92 

3712.66 

14 

EES 

9687.30 

mu 

KICKICI 

247.24 

0.0019 

0.9368 

1700.02 

0.0557 

9.92 

3712.08 

15 

EES 

9926.58 

ISIB 

249.49 

EESB 

0.9368 

1701.08 

0.0557 

9.92 

3711.64 

16 

BBI 

7078.00 

Ba 

EjB;! 

255.03 

0.0027 

0.9368 

1703.24 

0.0557 

9.91 

3710.76 

17 

EE13 

USB 

ktcWiil 

245.55 

0.0023 

0.9368 

1705.18 

0.0557 

3709.97 

18 

BBI 

3623.29 

tcBTil 

BtB 

282.53 

0.0035 

1714.97 

0.0555 

9.93 

3716.90 

19 

BBI 

EES 

3813.80 

BSEn 

285.65 

0.0034 

0.9421 

1716.30 

0.0555 

9.93 

3716.36 

20 

HUM 

BRU 

660.87 

liBl 

RIB 

403.98 

0.0287 

0.7954 

1425.80 

0.0596 

9.41 

3523.10 

21 

HBSl 

BBcl 

805.89 

BQ 

E!B 

367.91 

0.0219 

0.7954 

1433.35 

0.0594 

9.40 

3519.63 

22 

BBHI 

BBI 

792.95 

1.53 

40.35 

351 .90 

0.0185 

0.7954 

1430.43 

0.0595 

9.41 

3520.97 

23 

HES 

767.78 

kEISI 

323.00 

0.0174 

0.7954 

1428.73 

0.0595 

9.41 

3521.75 

24 

KIM 

BBI 

4254.54 

RBil 

346.65 

0.0041 

0.8307 

1495.67 

0.0585 

9.54 

3571.63 

25 

BEfl 

5041.56 

liltW 

RIB 

343.66 

0.0034 

0.8307 

1494.60 

0.0585 

9.54 

3572.11 

26 

n-g 

1869.69 

0.97 

42.40 

384.21 

0.0108 

0.8307 

1495.71 

0.0585 

9.54 

3571.62 

27 

3844.71 

BBI 

RBB 

370.15 

0.0051 

0.8307 

1501.07 

0.0584 

9.53 

3569.23 

28 

ttIB 

HBI 

1018.75 

liEltl 

BBI 

357.38 

0.0147 

0.8307 

1497.47 

0.0585 

9.54 

3570.83 

29 

1069.48 

BBI 

BBil 

460.05 

EEiai 

0.8484 

1536.23 

0.0579 

9.60 

3592.96 

30 

BB!!I 

Bra 

4844.67 

1.72 

39.04 

324.34 

0.0038 

0.8484 

1531.71 

0.0580 

9.60 

3594.94 

31 

ni«ti 

2320.81 

0.91 

37.15 

295.48 

0.0072 

0.8484 

1526.81 

0.0580 

9.61 

3597.09 

32 

nra 

1835.47 

1.83 

40.94 

373.83 

0.0091 

0.7335 

1371.76 

0.0605 

9.09 

3402.17 

33 

BBI 

3088.34 

BbSI 

345.43 

0.0050 

0.7335 

1375.28 

0.0604 

9.08 

3400.53 

34 

gfigl 

BB51 

3152.64 

332.71 

0.0047 

0.7335 

1373.63 

0.0604 

9.09 

3401.30 

35 

2670.25 

WM 

324.22 

0.0054 

0.7335 

1367.22 

0.0605 

9.09 

3404.30 

36 

^11 

BS\ 

2651 .76 

MM 

KHEEI 

323.04 

0.0054 

0.7194 

1341.25 

0.0610 

9.03 

3382.10 

37 

2779.43 

1^1 

298.74 

0.0048 

1343.85 

0.0609 

9.03 

3380.86 

38 

ssi\ 

KESII 

2133.71 

2.44 

37.03 

307.76 

0.0064 

0.7194 

1342  71 

0.0610 

9.03 

3381.40 

39 

EESII 

2054.57 

^91 

^^231 

283.57 

0.0062 

0.7194 

1344.60 

0.0609 

9.03 

3380.51 

40 

^1 

BSli 

1685.09 

248.21 

0.0066 

0.7194 

1348.95 

0.0609 

9.03 

3378.45 

Uest'-  estimated  particle  velocity  from  DELTAE,  (m/s) 
urms  -  root  mean  square  velocity,  (m/s) 

ReD  -  Reynolds  number  based  on  hydraulic  diameter 
LMTD  -  Log  Mean  Temperature  Difference,  (K) 
h  air  side  -  heat  transfer  coefficient  on  air  side,  (W/m^'K) 
j  air  side  -  j-Colbum  factor  on  the  air  side 

H20  vel  1  channel  -  velocity  of  water  in  one  heat  exchanger  channel 

Nusselt  -  Nusselt  number  for  water 

htubes  -  heat  transfer  coefficient  for  water  in  tubes 
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Appendix  III:  Table  of  Thermoacoustic  He-Ar  Properties 


0.000024 


0.000024 


EEB 

E13I 


Hot  Exchanger  Properties 


Thermoacoustic  Workina  Fluid  Prooerties 


Data  Point  rho  k  mu  DeltaK  DeltaN 


0.000063 


0.000063 


0.000024  0.000100  0.000063 


0.000024  0.000100 


0.000081  0.000052 


0.000024  0.000081  0.000052 


0.000024  0.000081  0.000052 


0.000025  0.000074  0.000049 


0.000025  0.000074 


0.000024  0.000091 


0.000047 


0.000025  0.000069  0.000047 


EE3inic!»niP»«WiMW3SSl| 

HBiliitjgiggBBSlSlHfSgMi 


iUItklHiHH 

i^gaEB^ _ 

It^giiESlkflfliaiaaiEhiiamiiffliBaa^EggEaiiBaiaiBi^^ 

iUJdlEEMll!il!!!!!!!H^BfiTATtV^iy6Mi^jffHii!iTrBlBBBgBn5fSSi! 

liu^EB^liniiimiMi^^^cgaigHiRHKgiajiB^^Illjilra 

itEggEEBgbiii!ai!iaEB55BifliiQii]iiaECTBMBl«N!«!t*i^4»^^^ 

ip4»iH!iBBBinin»«iini»Mi^«I«i»!«:i>4.-llffl!Bl«iasj|BBH^MMBBB^KK;?.-t;g 

in3aEBaaEB!«!«!«vapii?5MtK«i‘5EB[!ggi^^ 

inElBE3FlEEl»M«V-gH 
inE]BE5Bl 

|iimiciBB5{Rli»»»iwiw?li 

iKltilclCTiHHkmi»nnk-?li 
inWilclBBSH 


0.000101 


0.00002410.000101 


0.000101 


0.000101 


0.000101 


0.000101 


0.000101 


0.000101 


0.000024  0.000101 


||!l!!»!!!«!;aEggtB!cWtlBB!igg|E^ 

ItaiiiiaiiaEZfliiiiaaaia^^ 

IcaaMiaEgSuiiaaaiga^ 

IBBtri!i!SlBB3IB!t!iB1BB!ilclclBK!3ig 

icaaaiiiaiEEU!!^ 

litnnni«igiBBillB!3iBll«kit»!clciillriawM 


art  displace 


0.0027 


0.0022 


0.0012 


0.0035 


0.0022 


0.0031 


0.0035 


0.0044 


0.0044 


0.0040 


0.0041 


0.0038 


0.0043 


0.0052 


0.0053 


0.0046 


0.0052 


0.0034 


0.0035 


0.0004 


0.0005 


0.0007 


0.0007 


0.0023 


0.0028 


0.0010 


0.0021 


0.0007 


Hxler  Area 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


0.2307 


ReD 


1727.93 


Mj^BBBEl _ 

ifliiiaijEESIEI^iaigjli 

EiEg[!ii!iiijn=!-i*^»ii 

liB^BBSIHEHH;H 

nkniii»nmiiii3Ei!i»n];ii 

■ii«ii.tnBBI8HKw»i«i!l 

BiB!||^BIW^t!R>:jB 
Mj^BBBHI 
Wnf^BBBHlj 

iJIAt^CTBnSl 
mt.UlBBlgEl 
liBE3Br»MH 

■«ktt^BB!^i«k^J 

MtBn5ilB^EKE39 

lEEilEEBIBHEEEl 

laaisiflEEnim-liUEy 

iaaia!]ifl!iaflti'Iciik4--k 

ij^WiiiinnEnB?! 

i»I»{li{»lBB!in3i:ki.-tJJ 


1077.57 


2230.11 


679.41 


Hxref  Area  based  on  pd 


0.0680 


0.0553 


0.0313 


0.0871 


0.0556 


0.0777 


0.0872 


0.1086 


0.1085 


0.0999 


0.1024 


0.0938 


0.1067 


0.1283 


0.1319 


0.1133 


0.1282 


0.0836 


0.0883 


0.0113 
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0.0179 


0.0176 


0.0582 


0.0697 


0.0257 


0.0522 


0.0190 


0.0228 


0.0618 


0.0292 


0.0431 


0.0712 


0.0720 


0.0616 


0.0613 


0.0640 


0.0496 


0.0468 


0.0397 


Cold  Exchanger  Properties 


Data  Point 


B^EB^ 

EEEBSn 

BISlBBlBi 


mu 

DeltaK 

DeltaN 

0.000021 

0.000086 

0.000054 

0.000021 

0.000087 

0.000055 

0.000021 

0.000088 

0.000055 

0.000021 

0.000087 

0.000055 

0.000021 

0.000070 

0.000044 

0.000021 

0.000070 

0.000044 

0.000021 

0.000070 

0.000044 

0.000022 

0.000066 

0.000044 

0.000022 

0.000066 

0.000044 

0.000022 

0.000081 

0.000053 

0.000022 

0.000081 

0.000053 

0.000022 

0.000061 

0.000042 

HR3BB!!l3E^^i1E?r!;3 

HBlBlc!!BlBB!!53npSg 

FitElBSfUBBBBB 
PtBISSflBBlSIlBBSSi 


art  displace  1 

Hxfer  Area 

p 

Dh 

Dh  ReD  Hxref  Area  based  on  od 


BBIBilBBlilH _ 

BBtBilBBBBlBffiBni 

j^^BBBBlE!i!»ntiB 

BBflSlBBlBHBBgg 

BBBlilBBfiBIBiCTH 
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13 

[■lijcli!:] 

0.000022 

0.000061 

0.000041 

RBEI 

BISBBBI 

1.6667 

0.0054 

0.2307 

BISIBBl 

lifaMM 

7992.87 

0.1337 

14 

PfieliEl 

0.000061 

0.000041 

liEB 

BIBBBBI 

1.6667 

0.0065 

0.2307 

[iMMa] 

laxalail:] 

9687.30 

0.1610 

15 

0.000022 

0.000061 

0.000041 

EES 

^1X1 

BtiBBlH 

1.6667 

0.0067 

0.2307 

BKallifil 

BlalStlBI 

0.1653 

16 

BBBilSl 

0.000022 

0.000067 

0.000045 

BEB 

BliBBBI 

1.6667 

0.0057 

0.2307 

liXalljil 

Mililt:] 

0.1422 

17 

[il'lcMO 

0.000022 

0.000067 

0.000045 

EES 

[tkiltici:] 

1.6667 

0.0065 

0.2307 

Millslil 

ialaliit:] 

0.1613 

18 

BBISH 

0.000021 

0.000088 

0.000056 

BEBl 

nSBBlBI 

BlBBBBl 

1.6667 

0.0044 

0.2307 

ISEBIF:! 

3623.29 

0.1091 

19 

[BRl 

0.000021 

0.000088 

0.000056 

EES 

EliBBIi] 

BlBBBBl 

HBBS1 

0.0046 

0.2307 

BEngsi 

BEBIBI 

0.1151 

20 

EEliB 

tsxngsg 

0.000024 

0.000077 

0.000050 

EES 

iHiil 

liVAlc!! 

0.0006 

0.2307 

BlsnEl 

Mitatl:] 

0.0170 

21 

FEBl 

0.000023 

0.000075 

0.000049 

EES 

BIBBBig 

1.6667 

0.0007 

0.2307 

BKilTil 

MafiK:] 

805.89 

0.0196 

22 

tHHg 

0.000023 

0.000085 

0.000054 

BEBl 

ttilcMil 

BlBSBtS 

1.6667 

0.0009 

0.2307 

BlilBBl 

BfilitM 

792.95 

0.0233 

23 

i»E!gg:i 

0.000023 

0.000054 

BEBl 

KlltTiTl 

BEBBil 

1.6667 

0.0009 

0.2307 

talilMa] 

EPSBl 

767.78 

0.0230 

24 

0.000021 

0.000067 

0.000043 

EEO 

917.3 

BISBBEl 

0.0031 

0.2307 

Mallllil 

[•ialalia 

4254.54 

0.0766 

25 

RBElil 

liETOll 

0.000021 

0.000067 

0.000043 

EED 

917.3 

1.6667 

0.0037 

0.2307 

Mill-lil 

Laiai'ii.:] 

5041.56 

0.0910 

26 

BcBB 

BBEEI 

MIIIHHWI 

0.000067 

0.000043 

EES 

874.3 

BEBBBl 

1.6667 

0.0014 

laxalail:] 

0.0345 

27 

0.000022 

0.000066 

0.000042 

EES 

|:U|c| 

BliliBU 

1.6667 

0.0027 

0.2307 

MtlfAl 

lalalail:] 

3844.71 

0.0678 

28 

0.000023 

0.000077 

0.000049 

BEBl 

nSBBlBl 

BISBBS 

1.6667 

0.0010 

0.2307 

0.0160 

[axalaji:] 

1018.75 

0.0246 

29 

tBBl 

EE^l 

0.000022 

0.000083 

0.000052 

BEBl 

tBBBEl 

EIPP^ 

1.6667 

0.0012 

0.2307 

BBffiBl 

BfaBlBl 

EEIES 

0.0296 

30 

RHE51 

tiMlH 

0.000022 

0.000063 

0.000041 

BEBl 

817.4 

[ililtKl:] 

1.6667 

0.0032 

0.2307 

0.0160 

0.0018 

4844.67 

0.0798 

31 

0.000022 

MiV.MAl 

0.000041 

BEBl 

817.4 

[•lilijcia 

0.0015 

0.2307 

Mallji] 

laxalail:] 

BEO 

0.0391 

32 

CBH!1 

0.000021 

0.000088 

0.000056 

BEBl 

Ell¥3 

1.6667 

0.0023 

0.2307 

Milfilal 

Malail:] 

1835.47 

0.0565 

33 

nsRi 

fiEBB 

0.000021 

0.000088 

0.000055 

EES 

tBBtifil 

1.6667 

0.0037 

0.2307 

MaMal 

faWatt;! 

ESEES 

0.0928 

34 

BESE! 

0.000021 

0.000088 

0.000055 

BEBl 

BBWiM 

IBS 

0.0038 

0.2307 

MaEJtl 

EE3H 

3152.64 

0.0941 

35 

fHRI 

BBHH 

0.000021 

0.000088 

0.000055 

BEBl 

BBBBtil 

BEBBS 

1.6667 

0.2307 

HBIBBl 

BlaBlBl 

2670.25 

0.0805 

36 

BSBil 

0.000021 

0.000088 

0.000055 

BEBl 

HiBBW 

BfiBBS 

1.6667 

0.0032 

0.2307 

MaltAI 

|By 

2651.76 

0.0800 

37 

BEWM 

0.000021 

0.000088 

0.000055 

BEBl 

IBBBEl 

BEBBBl 

1.6667 

0.0033 

0.2307 

WaTBBl 

BBBHSI 

2779.43 

0.0834 

38 

[BE!!1 

BEH!3 

0.000021 

0.000088 

0.000055 

BEBl 

BBBBIil 

BEBBB1 

1.6667 

0.0026 

0.2307 

HBTBBl 

HBBIBI 

BBl 

0.0646 

39 

CHEEl 

BEEjjl 

0.000021 

0.000087 

0.000055 

BEBl 

IBBBBil 

BBBBBI 

1.6667 

0.0025 

0.2307 

HBIBBI 

HBBIBI 

2054.57 

0.0614 

40 

BEEjQ 

0.000021 

0.000088 

0.000056 

BEBl 

BBHtH 

BBBBBI 

1.6667 

0.0020 

0.2307 

lafaMH 

IBBSBH 

0.0514 

rho  -  density,  (kg/m^) 
k  -  heat  transfer  coefficient,  (W/m-K) 
mu  -  dynamic  viscosity,  (N-sW) 

DeltaK  -  thermal  penetration  depth,  (m) 

DeltaN  -  viscous  penetration  depth,  (m) 

Pr  -  Prandtl  number 

Cp  -  specific  heat,  (J/kg-K) 

beta  -  thermal  expansion  coefficient,  (K  *) 

gamma  -  ratio  of  specific  heats 

part  displace  -  particle  displacement,  (m) 

P  -  perimeter  of  heat  exchanger  fin,  (m) 

Dh  -  hydraulic  diameter  of  heat  exchanger,  (m) 
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Appendix  IV:  Table  of  Water-Side  Properties 


Hot  Exchanger  Properties 


water  properties 

Data  Point 

mu 

riSiRffl 

CP 

in 

■a 

1 

0.000926 

BRiSISSa 

4180.19 

n 

BEB 

2 

[iX01«l:l:il 

BEtl 

BER 

3 

0.000931 

IsltMTM 

ESESESl 

pni 

BEBl 

4 

0.000922 

0.000851 

miMH 

tini 

BEH 

5 

0.000927 

ithvm 

SBI 

BEB 

6 

0.000924 

0.000867 

vvimn 

4180.13 

BEn 

BER 

7 

0.000922 

nnPEi 

BBl 

BER 

8 

0.000927 

0.000874 

HR! 

BEB 

9 

0.000924 

0.000859 

liKMW'-l 

4180.13 

BTH 

BEH 

10 

0.000923 

0.000873 

4180.11 

EED 

BEB 

11 

MiIOIsWcl 

BUSES 

HRI 

BEB 

12 

BIBSEEEl 

HRI 

BEBl 

13 

BBiBsra 

railiWiil 

BEEXEl 

BIRI 

BER 

14 

HClBSra 

ntfiTO 

RRI 

BER 

15 

[iiimMc] 

BLliUd 

HRl 

BEB 

16 

nuloiiMci 

MtlililiWil 

PPVJUl 

4180.11 

BEH 

BEB 

17 

0.000921 

BFSES 

BEH 

BEB 

18 

0.000921 

[•11101:1.-1.1 

BEn 

BEB 

19 

r.Wiimi 

Enmi 

BEn 

BEB 

20 

0.000923 

PPlf/:U 

BKI 

BEn 

BER 

21 

0.000937 

MIIOMM 

mAtili 

BEfl 

BgR 

22 

[•nmotci 

RBPgg 

BEtl 

BER 

23 

mtmmn 

BEBEIH 

BEn 

BEB 

24 

rasiiliJjStn 

MHOIiliM 

EBl 

BER 

25 

PESra 

|iPVX:n 

Rpram 

SSI 

SS! 

26 

PEilSSSPl 

0.000915 

BESES 

BEn 

27 

HP555E!! 

KV/iM 

BEn 

ms 

26 

[Atniiiiwa 

0.000913 

RTligjl 

BEn 

BER 

29 

0.000933 

0.000910 

Eil:imc|.1 

BEn 

BEB 

30 

npitniiici.-i 

tiUtkoioi 

ni=T»gM 

BEn 

BER 

31 

iwi«ni!ic!;i 

0.000917 

k1:VA:!c1 

BESZS 

BEn 

BER 

32 

E!!VJ.-{| 

4179.51 

BEn 

BER 

33 

0.000813 

ESZSI 

4179.94 

BEn 

HBil 

34 

BEn 

Bni 

35 

l«A«IOI;Ml 

BEn 

BER 

36 

RK1R3 

BEn 

BER 

37 

EESSl 

4179.97 

BEn 

Bni 

38 

BEn 

BIBil 

39 

BEn 

BR»1 

40 

nkciiiiiii 

BQES 

BEn 

SB 

Cold  Exchanger  Properties 


Water  Properties 

Data  Point 

mu 

cp 

n 

Ei 

1 

0.000931 

0.000994 

ETHBH 

4180.32 

Bm 

sn 

2 

0.000933 

0.000994 

EBHB!! 

4180.35 

BEn 

S3 

3 

0.000935 

0.000994 

BiHBBI 

BEn 

S3 

4 

0.000929 

0.000994 

EBBKI 

BMS 

BEn 

BER 

5 

0.000933 

0.000994 

BBBI 

4180.37 

BEn 

S3 

6 

0.000932 

0.000994 

EBBBl 

4180.34 

BEn 

BEn 

7 

0.000931 

0.000994 

BiBBSl 

4180.32 

BEn 

S3 

8 

0.000935 

0.000994 

EEBRil 

4180.40 

BEn 

S3 

9 

0.000933 

0.000994 

EBBCT 

EILlilsS 

BEn 

BEB 

10 

0.000930 

0.000994 

4180.29 

BEn 

BER 

11 

0.000929 

0.000994 

4180.25 

BEn 

BER 

12 

0.000935 

0.000994 

4180.41 

BEn 

BER 

13 

0.000934 

0.000994 

4180.38 

BEn 

BER 
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14 

0.000933 

0.000994 

ESSBil 

4180.36 

BBn 

15 

0.000933 

0.000994 

4180.35 

cni 

EtH 

16 

0.000932 

0.000994 

4180.32 

uni 

EBl 

17 

0.000930 

0.000994 

4180.29 

HRI 

SIEI 

18 

0.000930 

0.000994 

4180.29 

HRI 

EfEl 

19 

Mmumti 

0.000994 

4180.27 

ESI 

ESI 

20 

0.000945 

0.000994 

4180.65 

HRI 

EISl 

21 

0.000940 

0.000994 

4180.53 

ESI 

EfH 

22 

0.000942 

0.000994 

4180.57 

ESI 

EIBI 

23 

0.000943 

0.000994 

4180.60 

ESI 

EiSEl 

24 

0.000941 

0.000994 

4180.55 

ESI 

EtEI 

25 

0.000941 

0.000994 

mil:!;] 

4180.57 

BBtl 

EtEl 

26 

0.000941 

0.000994 

4180.55 

ESI 

EB51 

27 

0.000937 

0.000994 

4180.46 

ESI 

EES 

28 

0.000940 

0.000994 

4180.52 

ESI 

SE3 

29 

0.000935 

0.000994 

4180.41 

ESI 

EEH 

30 

0.000938 

0.000994 

mil:!cl 

EUiliElil 

ESI 

SS 

31 

0.000941 

HUmiliMI 

4180.56 

ESI 

EER 

32 

0.000905 

0.000994 

4179.67 

PBn 

EIRil 

33 

0.000903 

0.000994 

4179.62 

ESI 

ESI 

34 

0.000904 

0.000994 

4179.64 

ESI 

EBil 

35 

0.000909 

0.000994 

BHBll 

4179.75 

Bfn 

EIBcl 

36 

0.000908 

0.000994 

4179.74 

HRI 

EISl 

37 

0.000994 

4179.70 

BBn 

Eni 

38 

0.000907 

0.000994 

4179.72 

BBn 

39 

0.000906 

0.000994 

4179.69 

BRI 

Eni 

40 

0.000903 

0.000994 

4179.61 

HRI 

EfEl 

mu  -  dynamic  viscosity,  (N-s/m^) 

mu(s)  -  dynamic  viscosity  at  surface  of  fluid,  (N-s/m^) 

density  -  density,  (kg/m^) 

k  -  heat  transfer  coefficient,  (W/m-K) 

Cp  -  specific  heat,  (J/kg-K) 

Pr  -  Prandtl  number 
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Appendix  V:  Table  of  Thermoacoustic  Raw  Data  for  the  MicroChannel  Heat  Exchangers 


Measurement 

Units 

Data  Points 

1 

2 

3 

pressure 

86 

l||^ 

KS 

Mix 

HE^ 

Freq 

Hz 

WliVAl 

Wac 

W 

Qc 

W 

Qh 

W 

Hrl'ii.ni 

IHHil 

DTh 

degC 

■EES 

DTc 

degC 

■EKS 

■IS 

Vrateh 

L/min 

I^SE] 

Vratec 

L/min 

KES 

stack  DT 

degC 

kPa 

12,20 

26.20 

18.80| 

Qc+Wac-Qh 

balance 

BD^ 

Thi 

degC 

EEESl 

Tci 

degC 

wtr-stk 

DT  -h  deg  C 

■EES 

wtr-stk 

DT  -c  deg  C 

piston  vel 

BBI 

■E^ 

pdr  phase 

BSIBI 

Tavph 

degC 

Tavgc 

degC 

MiMi\ 

lnbal/(Qc+Wac)*100 

MfJSi 

KH«H 

Appendix  VI:  Table  of  Thermoacoustic  Results  for  the  MicroChannel  Heat  Exchangers 


Hot  Exchanger  Results 


Data  Point 

Acoustic  Information  and  Calculation 

Water  Flow  in  HX  Tubes  Calculation 

nf?Bi 

rffuKi 

nami 

lOi 

h  air  side 

KiDiyiy 

H20  vel  1  channel 

friction  tactor 

I?™ 

1 

nBii 

HBCl 

530.27 

BEtl 

198.46 

0.0211 

0.33 

lcE!.i-gM 

0.04 

wmt 

2 

HRil 

1045.17 

jfcnif 

BKin 

151.19 

0.0081 

0.37 

KMtfikia 

0.04 

BH3RI 

3 

775.19 

EEa 

BiWJ 

212.79 

liBUgl 

0.16 

0.06 

UilikEU 

Cold  Exchanger  Properties 


IWflH 

1  ReD  1 

HTHti 

Bgl 

h  air  side 

1 1  air  side  1 

H20  vel  1  channel 

1  ReD  1 

friction  factor 

Nusselt 

htubes 

1 

Wcni 

iiUbJ 

BtlliM 

HBI 

39.81 

0.22 

0.05 

12.51 

1178.77 

2 

mi 

BBl 

BtIcM 

Eiig 

56.69 

0.19 

WWIMI 

0.05 

10.54 

50 


Appendix  VII:  Table  of  Thermoacoustic  He-Ar  Properties 


Hot  Exchanger  Results 


Data  Point 

IZ 

Thermoacoustic  Working  Fiuid  Properties 

liBl 

mu 

DeitaK 

DeitaN 

Pr 

EH 

ESSESI 

part  dispiace 

Hxfer  Area 

P 

Hxref  Area  based  on  pd 

1 

Hgl 

UmiriKI 

msm 

0.000152 

ImBBBftiSI 

1.6667 

0.0014 

0.3355 

DPra 

530.27 

0.0371 

2 

0.000112 

0.000071 

1.6667 

0.0015 

0.3355 

Flpra 

UiSklU 

3 

HBl 

0.000024 

0.000152 

0.000096 

1.6667 

0.0020 

0.3355 

vm 

775.19 

Cold  Exchanger  Properties 


liSRSijaiBi 

Thermoacoustic  Working  Fiuid  Properties 

liiBI 

mu 

DeitaK 

DeitaN 

Pr 

E3 

BfSfl 

aamma 

part  dispiace 

Hxfer  Area 

P 

ReD 

Hxref  Area  based  on  pd 

1 

0.000141 

ITBlil 

1.6667 

0.0018 

0.3355 

EBBiS 

817.26 

0.0486 

2 

EBEISlil 

0.000021 

8 

O 

o 

o 

d 

0.000063 

ITBlil 

1.6667 

0.0019 

0.3355 

0.0520 

Appendix  VIII:  Table  of  Water-Side  Properties 


Hot  Exchanger  Results 


Data  Point 

water  properties 

mu 

OffHl 

cp 

■9 

■a 

1 

bktuniiBiti 

4181.60 

BB51 

Efsn 

2 

0.000983 

ElhWAI 

PEil 

EHI 

3 

rai51!5E51 

Brnrsi 

EBB 

Cold  Exchanger  Properties 


Data  Point 

water  properties 

mu 

nwBi 

cp 

ra 

IS 

1 

ESI 

REnElil 

BB51 

EBB 

2 

0.000983 

BE13E1 

BbWAl 

BBlil 

^Qj, 
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Appendix  DC:  Thermoacoustic  Sample  Calculation 

This  sample  calculation  is  based  on  the  first  data  point  listed  in  the  first  four  appendices 
for  the  microchannel  rejection  heat  exchanger.  If  there  is  any  small  difference  caused  by  round¬ 
off  error,  the  correct  value  will  be  listed. 

•  Using  the  measured  data  from  Appendix  I,  construct  a  DeltaE  model  to  predict  the  particle 
velocity  in  each  heat  exchanger. 

■  M  =  2.9407  m/s 

•  Take  the  velocity  found  from  DeltaE  and  find  the  rms  velocity  by  dividing  by  yfl . 

U  2.9407  / 

•  Find  the  hydraulic  diameter 

■  Find  the  flow  area(Af.) 

-tL-  w,L, 

=  ;r  x0.0762^  -  0.0001016  x  12.272038  -  0.0018288  x  2.143235966 
=  0.01307508  m^ 

■  Find  the  total  heat  exchanger  surface  area(A) 

A  =  Af  +A,  =L,x2{w,  +t,)  +  Lx2(w  +  t) 

=  2.143235966  x  2(0.0018288  +  0.005969)  + 12.272038  x  2(0.0079375  +  0.0001016) 

=  0.230737332  m^ 

^  4LA,  4x0.0079375x0.01307508 

■  = - -  = - =  0.0018m 

A  0.230737332 

•  Find  the  Reynolds  number 

-  2.0794*11.029*.0018 

//  .000024 


52 


Find  the  log  mean  temperature  difference  as  the  difference  between  the  stack  end 
temperature  and  the  water  inlet  temperature. 


■  LMTD  = 


(T. 


W,0 


■T  ■) 

■*  W,l  > 


T  -T 

s  w, 

> - y 

T  -T 

I-*.!  ■‘W.oj 


.88 


T  n. 

-3.49 

[-2.61] 

=  3.03K 


•  Find  UA  by  dividing  the  hot  heat  exchanger  power  by  the  LMTD 


■  = 


Qh 

LMTD 


175.1 

3.03 


57.84 


•  Use  the  definition  of  UA  (Equation  6)  to  solve  for  ho-  Unfortunately,  the  hyperbolic  tangent 
in  the  efficiency  makes  it  impossible  to  explicitly  solve  for  ho.  Thus,  it  must  be  solved 
iteratively. 


UA  =  - 


■  +  Rh..  + 


tw 


hA 


1 


- ^ - +  6.387£:-5  + - ^ - 

.7799 *.231  3737.01  *.02646 


57.84 


where 


1973 

-  ri=\ - ^(l-;/.)  =  l - (1-. 7426)  =  .7799 

A  ^  .2307 

^  _  tanh(mL)  _  tanh(286.1*.0037)  _ 

~  ml  ~  286.1  *.0037 


.01608  */i„ 

190.4*8.065£-7 


=  286.1  m' 


■  R 


tw 


WDo/Dj) 

t-'tube 


=  6.387E-5  KAV 


with  the  iterative  solution  of 
■  =781.6  WW-K 
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